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I. INTRODUCTION 

With the precise measurements of the Z^ mass from LEP experiments, as well as other 
electroweak parameter measurements, the standard model of electromagnetic and weak in- 
teractions is being tested beyond tree level. Even without knowing the masses for the top 
quark and the Higgs boson, physical observables when compared to theoretical predictions 
constrain new physics. [|l| An improvement in the measurement of the W mass, currently 
with a value of 80.22 ± 0.26 GeV from CERN [I and Fermilab 0] Collider experiments, 
would even further those constraints. With the higher statistics expected in the 1993-1994 
Fermilab Tevatron Collider runs, the theoretical and systematic errors must be also be re- 
duced to make a significant improvement in the W mass measurement. Even though the 
dependence on the W transverse momentum (g^-) is not very strong in the W mass sample, 
the QCD predictions for the low qr production of IVs, including nonperturbative effects, 
must be better understood. 

The low qT behavior of VT's and Z's has been the subject of much theoretical work. 
The CoUins-Soper-Sterman formalism |^ for the low transverse momentum distribution in 
Drell-Yan has been applied to IV's and Z's by Davies, Stirling and Webber, Altarelli et 
al., [^ and the matching of low and high qx by Arnold and Kauffman. ^ A study of the 
parameterization of nonperturbative effects in the low qr region has recently been done by 
Ladinsky and Yuan. Theoretical uncertainties in the separate W and Z qT distributions 
are discussed in detail in Ref. P). Here, we discuss the relative uncertainties between the W 
and Z transverse momentum distributions. By considering the ratios 

and 

the theoretical errors are reduced because the uncertainties tend to cancel. By comparing 
these quantities with experimental results, one can test the theoretical treatment of low qx 



vector bosons. Alternatively, the ratios can be used in Monte Carlo models for the W mass 
measurement. 

The estimation of the theoretical uncertainties in R{qT, y) and R{qT) due to the parame- 
terization of nonperturbative effects is the main topic of this paper. In Section 2, we review 
the theoretical issues associated with the low transverse momentum behavior of W^s and 
Z's. Our numerical results are displayed in Section 3. Section 4 contains our conclusions. 
As all of the relevant formulae are collected in the Appendix of Ref . f^ , we do not include 
detailed formulae here. 

II. THEORETICAL ISSUES 

The calculation of the low transverse momentum behavior of Ws and Z's is not com- 
pletely straightforward because of the presence of large logarithms associated with powers 
of as- At low qt, the perturbative expansion in a^ is ruined because of the accompanying 
large logarithms. However, it is possible to resum the leading and next-to-leading ln{Q'^/q^) 
powers for each power of as to get an exponential factor. Collins, Soper and Sterman have 
done the resummation, f^ explicitly including transverse momentum conservation for mul- 
tiple parton emission at each order in a^. The momentum conservation appears through a 
Fourier transform involving the conjugate variable to qt, namely the impact parameter b. 
The resummation appears as a Sudakov-like exponent S, which is incorporated in a form 
factor, here called W. The quantity W can be computed perturbatively over the range of 
^QCD ^ 1/^ < Q- Schematically 

Wpertib, Q, Xa, Xb) ~ Yl ^^ij^) U I a{x A-, IJ,)fj/BixB, yU) exp[-S'pert (&, Q)] ■ (3) 

The parton distribution functions fi/A are labeled for parton i in hadron A. The full per- 
turbative form for W also includes convolutions of the parton distribution functions with 
other functions which are suppressed by as/{2TT). The coupling constants eij{V) depend on 
whether V = WotV = Z. 



Including a normalization parameter Ny^ V = W or Z, the differential cross section has 
the form 
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in the low qt limit. Eqn. 4 applies only to the low qx limit because not included is the part 
of the cross section which is regular as g^ ~^ 0. Arnold and Kauffman have shown that the 
regular part of the differential cross section for q^ and q^ has little effect below qr = 10 
GeV. As we shall see below, the nonperturbative contributions are sizable only for g^ ~ 10 
GeV, so we consider this region of low g^ and neglect the regular part throughout. The 
form factor W depends on the impact parameter, the mass of the weak gauge boson and 
on xa = {Ql \fS)ey and xb = {Q / \fS)e'y , where v^ denotes the hadron collider center 
of mass energy, and y is the weak gauge boson rapidity. The form factor also includes the 
dependence on the parton distribution functions, the resummed terms, a nonperturbative 
factor and the unphysical factorization scale /i. 

For h < l/Qi we roughly approximate W , since in the integral in eqn. 4, this range of 
h contributes little for qx ~ 10 GeV. However, for h ~ 1/ Kqcdi the perturbative expression 
for W breaks down. A cutoff in h is introduced in the standard way by 0,§|,H 

h 

b -^ K = 77^ (5) 
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where b^ax ~ VQo is characterized by a scale Qo at which perturbation theory is still valid. 
The substitution of 6 — > 6^ in Wpen is accompanied by a nonperturbative factor. Thus, W 
defined over the full range of b is written this way: 

W{b, Q, Xa, Xb, f-i) = t^pert(&*, Q, Xa, xb, /i) ■ e-^"''^'''^) . (6) 

By combining eqns. (3)-(6), the main theoretical uncertainties are evident: dependence on 
the parton distributions and Aqcd, the residual /i dependence, and the parameterization of 
the nonperturbative effects in Snp{b,Q). The nonperturbative effects dominate for g^ ~ 6 
GeV, so we focus our attention on the nonperturbative parameterizations. 



Not indicated in the formulae above is an additional dependence on three unphysical con- 
stants denoted Ci, C2 and C3 in Ref. [§]. Arnold and Kauffman showed that the dependence 
of da/dqT for W^s and Z's on these three parameters is weak for a range of parameters. |Q] 

The nonperturbative factor Snp has the form: ^ 

Snp{b, Q) = Gi{xA, Xb, b, bmax) + ^2(6, bmax) hl^Qbmax) ■ (7) 

In principle, Gi depends on the incoming parton type and is hadron dependent. Following 
Davies, Stirling and Webber , we will assume that Gi is independent of parton type and 
independent of xa and xb- Ladinsky and Yuan include a dependence on xa and xb in their 
parameterization, which we comment on below. The standard assumption is that Gi and 
G2 are proportional to 6^, P,|| and Snpip, Q) is written in terms of constants gi and §2'. 

Snp{b, Q) = [gi + 92 HQbmax/2))b^ . (8) 

This is, at large b, roughly equivalent to a Gaussian distribution of intrinsic momentum. 
Ignoring the b dependence in VFpert, one is left with the Fourier transform of a Gaussian, 
which is itself a Gaussian. The average transverse momentum squared in the Gaussian is 

{Qt) =4:[9l+g2HQbmax/2)]. 

Davies, Stirling and Webber §] have performed a fit to ISR (R209) and Fermilab (E288) 
fixed target data. Their numerical values for gi and (72 are: 

gi = 0.30 GeV^ c/2 = 0.16 GeV^ (npl) 
gi = 0.15 GeV^ ^2 = 0.41 GeV^ {np2) 
gi = 0.0 GeV^ ^2 = 0.60 GeV^ {np3) 

with bmax = 0.5 GeV^^. The nonperturbative parameterization np2 is their preferred fit. 
These parameters were used by Arnold and Kauffman in their analysis of theoretical errors 
in da/dqT- 

Evident from the discussion so far is the fact that the the distributions for Q = My/ and 
Q = Mz will not be very much different apart from the overall normalization factors and 



the combination of parton distribution functions, evaluated at slightly different values of xa 
and xb- We now proceed to consider the numerical results. 

III. NUMERICAL RESULTS 

In obtaining the numerical results, we have used M^ = 80 GeV and Mz = 91 GeV. The 
weak mixing angle was taken such that sin^ 9]y = 0.22. To compare nonperturbative param- 
eterizations, we have used the Harriman et al. |l^ parton distribution functions HMRSB 



with the four flavor MS scale Aqcd = 0.19 GeV, based on a fit to the BCDMS data. To 
compare the effect of a change in parton distribution functions on R{qT,y = 0), we also 
use the HMRSE distribution functions fit to the EMC data with Aqcd = 0.10 GeV. As in 
Ref. 0, the canonical choices of Ci = C3 = 26^""^ = bo and C2 = 1 are used. Here, je is 
the Euler constant. As indicated above, we neglect the regular part of the differential cross 
sections, which contribute on the order of less than 1% in the range of qt < 10 GeV. §] All 
of the figures show results for the Tevatron Collider at viS = 1.8 TeV. Similar results obtain 
for the CERN Colhder. 

A. Nonperturbative parameters gi and (72 

The largest uncertainty for qx ^ Q GeV comes from the choice of parameters gi and g2 in 
eqn. (8). Figure 1 shows the distribution da / (dqTdy) at y = as a function of qx of the W, 
including both charges of the W in the rate. The three predictions vary from the average 
of the three curves by as much as 40% for very low g^. For g^ ~ 6 GeV, the the sensitivity 
to the nonperturbative parameterization is significantly less. Unfortunately, it is precisely 
below ~ 6 GeV where the differential cross section is the largest. A similar uncertainty is 
obtained for the Z^ transverse momentum distribution as a function of the nonperturbative 
parameterization. 

Because the W and Z masses are close and there is only log dependence on the masses, 
the nonperturbative parameterization Snp{b,Q) and perturbative exponent Spej.t{bi,,Q) are 



nearly identical for Q = Mw and Q = Mz- One expects that the ratio of the W to Z 
transverse momentum distributions will be approximately constant for small qt- Indeed, 
this is the case. Figure 2 shows the relative size of the W to Z distributions aX y = for 
npl — np3. Qualitatively, the size of R{qT, 1/ = 0) can be understood from a combination of 
coupling constants and parton distribution functions. Consider the quantity 

Dv{xA, xb, /x) = ^ el{V)fi/p{xA, ^)fj/p{xB, /i) (9) 

By taking 

Dw{Mw/^, Mw/^, Mw)/DziMz/VS, Mz/VS, Mz) , (10) 

one obtains the value of 3.04. If one were to evaluate the W and Z parton luminosities 
at the same values of xa and xb, the ratio would equal 2.34, as a consequence of their 
different coupling constants. Deviations of the npl — np3 curves from the average value 
Ra.vg{(lT,y = 0) are less than a percent. 

To understand the increase in R as qt decreases, it is useful to look at W for W^s and Z's, 
which is independent of qt- In Figure 3a, we show the form factors W{b,Q,XA,XB,bo/bi^) 
with Xa = Xb = Q/yS for Q = M\y and Q = Mz, and with the nonperturbative parameters 
np2, as a function of ln(6 ■ GeV). The low b behavior is an extrapolation of the value of 
W aX b = 1/Q, to lower values of b. The integral over b oi W also includes a factor 
of bJolbqT), so the differential distribution is fairly insensitive to such small values of b 
where the extrapolation must be performed. For large values of b, the nonperturbative 
parameterization cuts off the integral. Figure 36 shows the ratio of the form factors W, for 
the three sets of nonperturbative parameters. We cut off the range of ln(6 ■ GeV) displayed 
because at low b, the shape of the ratio is an artifact of the low b extrapolation. The curves 
for npl — np3 separate at 6 ~ lO^^GeV^^, because of the mass dependence in Snp, with the 
largest ratio belonging to np3. 

At low qt, the average value of b is for "large" b. For example, for W^s at qx = 0.5 
GeV, (b) = 0.8 GeV^^. As qr increases, the Bessel function in the integral oscillates more 
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rapidly as a function of b, and (6) decreases, moving into the perturbative range. Thus, as qr 
increases, the sensitivity to nonperturbative physics decreases and the calculation is more 
reliable. Figure 36 also shows the trend that -R(q'r) and R{qT, y) decline with increasing 
qt- If one trusts the form of the nonperturbative parameterization, then the uncertainty in 
R{(lT,y = 0) is less than a few percent. Even if one estimates the uncertainty in R by the 
deviation from the ratio of coupling constants and parton distribution functions, one is still 
left with a small uncertainty due to nonperturbative parameterizations of order 5%. 

Figure 4 shows the ratio Riqr) (where now the rapidity of the gauge boson has been 
integrated) for npl — np3, as a function of qx- The shapes of the distributions are nearly 
identical to those of Figure 2, however, with an increase in the normalization. This is due 
to the fact that the ratio of the W to Z parton luminosities increases with increasing \y\. 
We have the same qualitative conclusions about the uncertainty in the ratio Riqr) as in the 
y = case. 

B. The Ladinsky-Yuan parameterization 

The Ladinsky-Yuan parameterization [§ is characterized by three parameters and de- 
pends on S as well as Q. In a recent reanalysis of the R209 and E288 low qx Drell-Yan data, 
Ladinsky and Yuan (LY) have postulated a modified form for Snp'- 

Snp{b,Q,XA,XB) = gib^ + gig3blia{100xAXB) + g2b^i^{Qbmax/l-G) (H) 

The inclusion of the xa and xb dependence is to account for changing (gf,) as a function of 
Q'^/S. The W and Z qx distributions have peaks at slightly higher g-p values than in the 
two parameter cases. We use their central values of (/i, (72 and g-^ to characterize the effects 
of this parameterization on R{qT,y) and R{qT)- Their central values are gi = 0.11 GeV^, 
g2 = 0.58 GeV^ and g^ = —1.5 GeV^^, and their fit was performed using the CTEQ parton 



distribution functions. |Tl]] As we are interested primarily in the ratios, we continue to use 



the HMRSB parton distribution functions. The LY curve for R{qT) at v^ = 1.8 TeV is 



shown as the dotted hne in Figure 4. The deviation from the two parameter result is less 
than a few percent. 

C. Parton distribution dependence 

Figure 5 shows the ratio of the W and Z qt distributions for two different parton distri- 
bution functions. We show R{qT,y = 0) for the HMRSB and HMRSE parton distribution 
functions with np2. The HMRSE distributions yield the slightly lower value of R, but the 
two predictions for i?avg are within 2% of each other at fixed qt- 

IV. CONCLUSIONS 

The values of gi and g2 determined by Davies et al. are clearly not the best cur- 
rent values because of the advances in our knowledge of the parton distribution functions. 
Ladinsky and Yuan, [^ by refitting the data with current parton distributions, have made 
a significantly better choice. Nevertheless, uncertainties still exist in the form of the pa- 
rameterization, affecting the extrapolation from Q ~ 5 GeV to Q ~ 80 — 91 GeV. Parton 
distribution uncertainties at large Q compound the uncertainties. |^ 

We have shown that by taking the ratio oi W to Z qt distributions, the theoretical 
uncertainties due to the parameterization of nonperturbative effects tend to cancel. For the 
range of qx < 10 GeV, the ratio R^Qt) is ~ 3.3, slightly higher or lower depending on qx- 
When one takes the ratio of the total cross sections obtained by integrating the resummed, 
matched differential cross sections for np2, one obtains R = 3.27. R is less than 3.3 
because with increasing qt, -R(q't) declines roughly linearly, with R{qT = 100 GeV) = 2.66. 
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In our focus on the nonperturbative effects, we have neglected the regular part of the 
differential cross section which comes in at the percent level. A full study would also include 
an analysis of the other parameter dependence described in the text, e.g., Ci, C2, C3, bmax, 



etc. We expect, however, that in the ratios R{qT, y) and Ri^qr) the dependence should largely 
cancel. 
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FIGURES 
FIG. 1. The distribution da/{dqT<iy) at y = as a function of qx of the W for the sum of IF"*" 
and W' production at v5^ = 1.8 TeV, with nonperturbative paranieterizations npl (sohd hne), 
np2 (dashed line) and np3 (dot-dashed hne). 

FIG. 2. R{qT, y = 0) defined in eqn. 1 for npl — np3, as a function of qx for \/S = 1.8 TeV. 

FIG. 3. a) IFpert(^*, Q) exp(— S'np(6, Q)) with y = 0, vS = 1.8 TeV and fi = bo/K, as a function 
of ln(6GeV). The sohd hne is for Q = M\y and the dashed hne is for Q = Mz- b) Ratio of W for 
W^ to W for Z^ production, as a function of ln(6GeV), for npl (sohd hne), np2 (dashed line) and 
np3 (dot-dashed hne). 

FIG. 4. R{qT) defined in eqn. 2 for npl — np3, (sohd, dashed and dot-dashed hues) as a function 
of qx for v^ = 1.8 TeV. The dotted hne indicates Riqr) for the Ladinsky-Yuan parameterization. 



FIG. 5. R{qT, y = 0) for the HMRSB (sohd hne) and HMRSE (dashed hne) parton distribution 
functions for V^ =1.8 TeV and np2. 
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